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ABSTRACT 

G-band bright points (GBPs) are thought to be the foot-points of magnetic 
flux tubes. The aim of this paper is to investigate the relation between the dif¬ 
fusion regimes of GBPs and the associated longitudinal magnetic held strengths. 
Two high resolution observations of different magnetized environments were ac¬ 
quired with Hinode /Solar Optical Telescope. Each observation was recorded 
simultaneously with G-band hltergrams and Narrow-band Filter Imager (NFI) 
Stokes I and V images. GBPs are identihed and tracked automatically, and then 
categorized into several groups by their longitudinal magnetic held strengths, 
which are extracted from the calibrated NFI magnetograms using a point-by- 
point method. The Lagrangian approach and the distribution of dihusion indices 
(DDI) approach are adopted separately to explore the dihusion regime of GBPs 
for each group. It is found that the values of dihusion index and dihusion co¬ 
efficient both decrease exponentially with the increasing longitudinal magnetic 
held strengths whichever approach is used. The empirical formulas deduced from 
the htting equations are proposed to describe these relations. Stronger elements 
tend to dihuse more slowly than weak elements, independently of the magnetic 
hux of the surrounding medium. This may be because the magnetic energy of 
stronger elements is not negligible compared with the kinetic energy of the gas 
and therefore the hows cannot perturb them so easily. 

Subject headings: techniques: image processing — Sun: photosphere — methods: 
data analysis — methods: statistical 
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INTRODUCTION 


Magnetic fields are ubiquitous in the photosphere and interact with convective 
plasmas. Photospheric convective plasma flows advect magnetic helds, and eventually 
concentrate them. This interaction may cause magnetic reconnections and the excitation 


i Alfven 

1947; 

Parker 

1998 


r uehes et al. 

2993 

; De Pontieu et al. 

2997; 

Tomczv 

c et al. 

2997; 

van Balleeooiien et al 

2911 

Stanealini et al. 

2913a 

b; 

Giannattasio et al. 

2914 

). Thus, it is 


impor tant to investigate the interaction betwee n small-scale magnetic elements and plasma 


flows flHughes et ah 


2003 


Viticchie et ah 


20061) . Given that magnetic helds are passively 


transported by plasma hows, the motion of magnetic elements is a manifestation of intrinsic 
organizing processes and can be described in terms of dihusion. 

Bright points in the photosphere are thought to be the foot-points of the magnetic 


hux tubes tha' 


( Stenho 


1985 


the convecti ve motions of granules push violently into intergranular lanes 


Splanki 


19931) . G-band observations are suitable for studying bright points 


because they appear b righter due to a re duced abundance of the GH-lines molecule at 


a higher temperature fjSteiner et al 


200l|), and are often referred to as G-band bright 


points (hereafter GBPs). Therefore, it is possible to measure the dynamics of photospheric 
magnetic hux tubes, although only a fraction of magnetic elements are thought to be 


associated with GBPs bv observations, maene 

:o-hv( 

i rodvnamical simulations, or semi- 

empirical models o 

hux concentrations ( 

Keller 

1992; 

Berger & Title 

2991; 

Sanchez Almeida 

o 

o 

C^l 

Steiner et al. 

2991; 

Schiissler et al. 

2993; 

Garb 

Bson et al. 

2994 

Shelvag et al. 

2994 



Beck et al. 


29971 : 


Ishikawa et al. 


2997 


de Wijn et al. 


29981). 


In the past decades, many authors focused on the photospheric dispersal of magnetic 
elements or bright points in a held of view (FOV), such as an active region (AR) or a 
quiet Sun (QS) region. Dihusion processes represent the efficiency of dispersal in the 
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photosphere, which uses a diffusion index, 7 , quantifying the transport process with respect 
to a normal diffusion (random walk). Historically, normal-diffusion (7 = 1 ) is the hrst 
known diffusion process. It characterizes a trajectory that co nsists of su c cessive random 
steps and is descri 


Lemons fc Gvthiel 


jed b y a simples t forni of diffusion theory (iFick 


1855 


Einstein 


1905 


199711 


LeightonI (1196411 hr st held that the dispers al rat e of magnetic 


region s in the photosphere is normal. Later, 


Jokipii fc Parker! (1196811 and 


Muller et ah 


(1199411 suggested that photospheric magnetic elements move in a random walk due to 
supergranular hows. When 77 ^ 1 , the process is termed anomalous dihusion. The motion 
of magnetic elements in ARs or network regions in the QS is sub-dihusive (7 < 1) due to 


being trapped at stagnation points (i.e. points with nearlv 

0 zero horizontal velocitv: sinks 

of how held: Lawrence & Schriiver 

1993 

Simon et al. 1995: 

Cadavid et al 

199c 

)). More 

evidence based on hig 

super-dihusive (o' > 1 ; 

1 resolution 

at a s 

low that the photospheric BP motion is probablv 

Berger et al. 

1998: 

Lawrence et al. 2001: Kevs et al. 

2014; 

Yang et al. 


2015 


Typically, 


Abrameuko et al 


(1201111 indicated that the 7 value increases from a plage 


area in AR (7 = 1.48) to a QS region (7 = 1.53), and to a coronal hole (CH; 7 = 1.67). 


The dihusion coefficient, K, expresses the rate of increase in the dispersal area in 


units of time for magnetic elements or GBPs. The iL values range from 60 t o 176km^s ^ 


in A : 

Is or network regions in the OS (Wang 1 

988; S 

chriiver & Martin 

1990: Berger et al. 

1998; 

Hagenaar et al. 

1999 


Biannattasio et a 

. 2014 

; Kevs et al. 

2014) 

, and from 190 to 

400 km^ s ^ in QS regions 

S 

chriiver & Martin 

1990; 

Berger et al. 

1998 

: Giannatt 

asio e 

t al. 

2014; 

Jafarzadeh et al 

2014: 

Yang et al. 2015 

). In particular. Abramenko et al. 

(2011 

) 


compared the K values in an AR and QS region, which were about 12 and 19km^s 
respectively, whereas it was 22 km^ s“^ in a CH region. 


However, studies focusing on the dihusion of GBPs at diherent magnetic held strengths 
are scarce. Recent high spatial and high temporal magnetograms acquired with the Hinode 
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/Sola r Optical Telescope (SOT: 


Kosuei et al. 

— 

o 

o 

CVl 

Ichimoto et al. 

2008; 

Suematsu et al. 


20081 ) provide an unpr ecedented oppor 


strength. For instance, 


Utz et al. 


:unity to map a GBP to its co-spatial magnetic field 


(120131) extracted the corresponding magnetic field strength 


of each GBP from the longitudinal magnetogram within the spectro-polarimeter (SP) data. 
Nevertheless, high temporal SP data have such a narrow FOV that they are unsuitable for 
tracking the complete trajectories of GBPs. Instead, Stokes I and V images of Narrow-band 
Filter Imager (NFI) are co-spatial and co-temporal with G-band images, and keep a large 
FOV in observations. Thus, it is feasible to extract the simultaneous longitudinal magnetic 
field strength of each GBP during its lifetime. The aim of this paper is to study the relation 
between the dispersal of GBPs and the associated longitudinal magnetic field strengths. It 
will shed light on how magnetic elements of different longitudinal fields diffuse on the solar 
surface, and then assist the study of interaction between convections and magnetic fields. 


The layout of the paper is as follows. The observations and data sets are described 
in Section 2. The data reduction and analysis are detailed in Section 3. In Section 4, the 
relation of diffusions of GBPs on different longitudinal magnetic field strengths is presented. 
Finally, the discussion and conclusion are given in Section 5 and 6, respectively. 


2. DATA SETS 

We used two data sets of different magnetized environments acquired with Hinode 
/SOT. Each data set comprises G-band filtergrams (BFI) and Stokes / and V images (NFI). 
The G-band data with a wavelength of 4305 A are suitable for photospheric bright point 
sensitive investigations. The circular polarization I and V images can measure longitudinal 
magnetic fields. 

In data set I, there is an AR (NOAA 10960) between 2007 May 30 and June 14. 
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The region produced numerous C- and M-class flares during that period. We adopted a 
timeseries of G-band images taken between 22:31:21 and 23:37:26 UT on June 9. The images 
were obtained at a 30 s cadence with an exposure time of 0.15 s. The spatial sampling 
is 0.109" over a 111.6"x 111.6" FOV. The center of FOV is a; = 429.5" and y = -158.5” 
corresponding to a heliocentric angle of 28.2° equaling a cosine value of 0.88. We also used 
a co-spatial and co-temporal timeseries of Stokes / and V images taken 200 mA blueward of 
the Na l D 5896A spectral line with a spatial sampling of 0.16" over a 327.7"x 163.8" FOV. 

Data set II was recorded on 2007 March 31. It covers a quiet Sun region near the 
solar disc center. The FOV of the G-band data is 111.6"x55.8" with a resolution of 0.109", 
and its center pointed to solar coordinates of a; = 207.7" and ?/ = -130.3". This position 
corresponds to a heliocentric angle of 12.4° equaling a cosine value of 0.98. The data set 
starts from 11:36 UT until 12:40 UT, with a temporal sampling of 35 s. The timeseries of 
co-spatial and co-temporal Stokes I and V images taken 120 mA blueward of Fe l 6302.5A 
spectral line has an FOV of 163.8"x81.9" with a resolution of 0.16". 


3. DATA REDUCTION AND ANALYSIS 

The flltergrams were calibrated and reduced to level-1 using a standard data reduction 
algorithm fg_prep.pro. The projection effects of both data sets were corrected according to 
the heliocentric longitude and latitude of each pixel. 


3.1. Alignment 

The G-band images and the Stokes I and V images was aligned carefully. The 
temporally closest images were chosen first, and the different spatial sampling were 
overcome by bicubic interpolating the Stokes I and V images to the corresponding spatial 
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sampling of the G-band images. A fter that, a sub_pixel level image registration procedure 


(Fene: et ah 

2QI2; 

Yansf et ah 

2015) 


201511 was employed for spatial alignment; All G-band images 


in the timeseries were aligned to its hrst image, and then the Stokes I and V images were 
aligned to the G-band images according to the displacement between the simultaneous 
Stokes I image and the G-band image. Finally, all of the images were cut for the same 
FOV. 


3.2. Calibration of NFI Magnetogram 


The NFI Stokes I and V images of two data sets were measured with Na l D 5896A and 
Fe I 6302.5A spectral lines, respectively. Galibration of the NFI magnetograms is needed 
because they do not allow a full Stokes inversion to be performed. We calibrated the NFI 
longitudinal magnetic held strengths with a reference to the longit udinal magnetog rams 


in the level-2 SP data, which is inverted from the MERLIN code (ILites et ah 


20071). As 


indicated in the previous studies for magnetic helds outside sunspots, a linear relation. 



The NFI V/I images were desampled to consist with the spatial sampling of the SP 
longitudinal magnetogram. Then the narrow stripes of the timeseries NFI V/I images 
were cut and aligned with a reference to the co-spatial and co-temporal stripes of the SP 
longitudinal magnetograms by image registration techniques. The SP observation of data 
set II was simultaneous with the NFI images of Fe spectral line, so the calibrated coefficient 
of data set II was measured using the simultaneously SP data directly. Unfortunately, 
co-spatial and co-temporal SP data were not taken with the NFI images of data set 
I, so another NFI observation (the same as Na spectral line of data set I, which had 

























simultaneously SP data), was employed to measure the [3 value of Na line. This observation 
was recorded from an AR between 20:20:05 and 20:59:29 UT on 2007 July 1, which was 
temporally close to data set I. Both SP data were obtained by scanning a narrow slit over 
an area of 297.1"xl63.8" with a resolution of 0.30"x0.32" per pixel. 

Figure [T] shows the plots of NFI V/1 versus Bgp of Na and Fe spectral lines, respectively. 
The SP longitudinal magnetograms were multiplied by corresponding hlling factors. The 
pixels where the values were less than noise of SP longitudinal magnetograms or NFI V/I 
images were excluded. Only the pixels whose absolute field strengths are less than 1000 G 
were taken into account in data set I. With linear regression analyses of the relations of Na 
and Fe spectral lines, the (3 values are 13.1 and 35.6 kG with correlation coefficients of 0.95 
and 0.87, respectively. Gonsequently, longitudinal magnetic fields at any location in the NFI 
data of data set I and II were determined by multiplying the V/I with the corresponding (3 
values. The mean longitudinal magnetic held strengths of data set I (excluding the sunspots 
region) and II are 132 and 64 G, respectively. In addition, we selected a relatively quiet 
region (about l"xl") in the NFI magnetograms to quantify the noise of NFI magnetograms. 
The standard deviations (cr) are 10 and 20 G, respectively. 

For estimating the variation of the calibrated coefficients in different observations, 
we also processed some other data sets observed with Na l D 5896A lines. The result is 
that their (3 values are stable and limited in the range of 10%. Therefore, it is feasible to 
calibrate the NFI magnetograms of data set I using the calibration coefficient of the Na 
spectral line measured by another observation. 
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3.3. Detecting and Tracking GBPs 


A Laplacian and morphological dilation algorithm flFeng et ah 


201311 was ns ed to detect 


GBP s in each G-band image, and then a three-dimensional segment algorithm flYang et ah 


201411 was employed to track the evolntion of GBPs in the timeseries of G-band images. 


Fignre |2] shows a G-band image of data set I and the corresponding NFI magnetogram, 
in which the positions of the GBPs are highlighted in red. The GBPs cover 3.5% of the 
selected FOV. Fignre IS show s the images of data set II, bnt the GBPs only cover 1%. 


Sanchez Almeida et ah 


regions. Later, 


(120041) proposed tha 


Sanchez Almeida et al. 


the GBPs cover 0.7% in the internetwork 


(120101) detailed that the valne was between 0.9% and 
2 .2% in QS regions. Recently, the fractional a rea occnpied by GBP s in ARs is measnred as 


twice to triple larger than that in QS regions (iRomano et al. 


2012 


Feng et al. 


201311 . 


For redncing detection error, these GBPs were discarded if (1) their eqnivalent 
diameters are less than 100 km or greater than 500 km, (2) their lifetimes are shorter than 
60s, (3) their horizontal velocities exceed 7kms“^, (4) their lifecycles are not complete, or 
(5) they merge or split dnring their lifetimes. As a resnlt, a total of 103,023 bright points 
remain in 132 images, yielding 18,494 evolving GBPs of data set I; 19,349 bright points 
remain in 110 images, yielding 3,242 evolving GBPs in data set II. 


3.4. Extraction of the Magnetic field Strengths of GBPs 


The peak longitndinal magnetic held strength of each bright point was extracted from 
the corresponding NFI magnetogram after its region was identihed. The reason why a peak 
valne was adopted rather than the mean or median is that the scales of GBPs are so small 
that slight errors in image alignment c an degrade average or median valnes signihcantly 


(iBerger fc Ronppe van der VoortJ 120071 ). 
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Fig. 1. — Plots of NFI Vjl versus of different spectral lines. Panel (a): plot of NFI Vjl 
versus of 200 mA blueward of the Na l D 5896A spectral line. The data were recorded 
between 20:20:05 and 20:59:29 UT on 2007 July 1. Panel (b): plot of NFI V/I versus Bgp of 
data set II, with 120 mA blueward of the Fe I 6302.5A spectral line spectral line. See the 
electronic edition of the Journal for a color version of this hgure. 
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Fig. 2. — Panel (a): a G-band image of data set I (NOAA 10960), which the size is 
104.l"x 105.9". Panel (b): the co-aligned NFI magnetogram; the GBPs identihed from the 
G-band image are highlighted with red. The x and y coordinates are ticked by arcseconds. 
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Fig. 3.— Panel (a): a G-band image of data set II (a quiet Sun region), which the size is 
109.8"X54.3". Panel (b): the co-aligned NFI magnetogram; the GBPs identified from the 
G-band image are highlighted with red. The x and y coordinates are ticked by arcseconds. 
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The absolute strongest magnetic field strength, during the evolution of each 
GBP was calculated because it represents the peak state. Taking the noise of the NFI 
magnetogram into account, we discarded the GBPs from data set I and II with a B are 
below 50 and 100 G (about 4 ct), respectively. About 97% GBPs of data set I fall into a 
range of 50-1000 G, while 97.5% of data set II fall into 100-450 G. Gonsequently, these 
GBPs of data set I and II were categorized into 19 and 7 groups by setting the same bin as 
50 G, respectively. 


3.5. Diffusion Approaches 


A traditional analysis of the diffusion process is the Lagrangian approach, which is 
efficient when we deal with diffusive properties of tracers within turbulent fluid flows 


f|Monin &: laglom 


19751) . The approach consists of two main steps: (1) computing the 


spatial displacement, AZ, of an individual GBP as a function of time interval, r, measured 
from its first appearance; ( 2 ) calculating the mean-square displacement, ((A/)^), for each 
time interval as a function of r. The power index, 7 , of the mean-square displacement is 
defined as 

{{Alf) = Ct\ (1) 

where, C is the coefficient of proportionality. Usually, 7 and C are derived as the slope and 
the exponential of the intercept on the y axis of the spectrum over a range of r on a log-log 


scale, r espectively. However, this approac 


GBPs flDvbiec &: Gudowska-Nowak 


2009; 


1 is not ideal for study ing the random motions of 


Jafarzadeh et af 


20141) . These authors indicated 


that the second step might cause the mixing of different diffusive processes. Thus, they 
suggested that the square displacement of each GBP ought to be calculated separately. The 
7 value could be measured by the slope of its own square displacement on a log-log scale. 
Then, the mean diffusion index and the associated standard deviation (the square root of 
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the variance) could be obtained by fitting the distribution of the 7 values. We named this 
approach distribution of diffusion indices (DDI). In this study, both of the approaches were 
adopted and compared. 


4. RESULT 


4.1. Lagrangian Approach 


The mean-square displacements of GBPs in different longitudinal magnetic held 
strength bins versus the times of data set I and II are displayed on a log-log scale in Figure H] 
(a) and (b), respectively. Many previous authors have proposed that the estimation of 
7 using the Lagrangian approach may be strongly affected by timescales. In particular, 


Abramenko et ah 


(120111 1 interpreted that the 7 value measured for small timescales 


< 


300 s) could represent the intrinsic diffusion of GBPs. As shown in Figure HI the lengths of 
mean-square displacements are different and the tails are not regular because only a few 
GBPs have long lifetimes. If the 7 value was calculated by htting the whole mean-square 
displacement, it would be determined by a few long-lived GBPs. Therefore, we analyzed the 
mean-square displacements for the timescale r < 300 s to probe the intrinsic depiction of 


the diffusion index. The 7 values and the associated standard deviations were conveniently 
captured, respectively, by the slopes of mean-square displacements for the timescale r < 
300 s on a log-log scale within 95% conhdence intervals. 


In Figure El (a), the mean diffusion indices and the associated standard deviations of 
the different longitudinal magnetic field strength bins of data set I are illustrated using error 
bar in black solid line. Taking the quantity of each bin as a weight, we analyzed the relation 
between the diffusion indices and the longitudinal magnetic held strengths using a weighted 
curve htting. It hts well with an exponential function (black dashed line). The 7 value 
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Fig. 4.— The mean-square displacement ((A/)^) of GBPs as a function of time, r, on a log- 
log scale in different longitudinal magnetic field strength bins by the Lagrangian approach. 
Panel (a): The mean-square displacement of data set I, which the longitudinal magnetic held 
strength bins range from 50 to 1000 G. Panel (b); The mean-square displacement of data 
set II, which the longitudinal magnetic held strength bins range from 100 to 450 G. 
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and the standard deviation are 1.61±0.17 (snper-diffusion) inside the bin of 50-100G. 

As the longitudinal magnetic held strength increases, the 7 value decreases. At a strong 
longitudinal magnetic held strength, the gradient becomes small and the 7 value arrives at 
~1.00 (normal-dihusion). The empirical formula deduced from the htting equation is given 
to be: 

7 (S) = + c, ( 2 ) 

where B is the longitudinal magnetic held strength in kG. The parameters a, b, and c are 
0.77±0.10, -1.95±0.91, and 0.96±0.15 under the 95% conhdence interval, respectively. 

The dihusion indices of data set II are also shown in Figure |5] (a) with a red dotted 
line. Gompared with data set I, the trend is simple because the longitudinal magnetic held 
strengths of GBPs only range from 100 to 450 G. The 7 value continues to mostly decrease 
from 1.70±0.07 to 1.27±0.10, except the hrst 7 value is smaller than the second, and the 
hfth is slightly smaller than the sixth. In Figure 0] (b), it can be seen that a sudden drop of 
the mean-square displacement happens at r ~ 300 s in the hrst longitudinal held strength 
bin. This leads to a small 7 value with a large standard deviation. Because of the limited 
range of longitudinal magnetic held strengths, we neglected the htting formula. 


We then established 

;he di 

formula ( 

Monin & laelom 

1975) 


ihusion coefficient, K, of anomalous dihusion with the 




K{t) = —T 


( 3 ) 


where 7 , r, and C are deduced from equation ([T]). Figure [5] (b) shows the relations between 
K and B for the timescale r < 300 s. The K value of data set I decreases exponentially 
from 143±50 to 26±4km^s“^. The empirical formula is given as: 


k{B) = + c. 


( 4 ) 


where B is the longitudinal magnetic held strength in kG. The parameters a, b, and c are 
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Fig. 5.— Panel (a): the relations between the diffusion indices of GBPs and the longitudinal 
magnetic held strengths for the timescale is less than 300 s by the Lagrangian approach. 
Panel (b): the relations between the diffusion coefficients of GBPs and the longitudinal 
magnetic held strengths for the timescale is less than 300 s by the Lagrangian approach. The 
mean values and the standard deviations of data set I are illustrated using the error bars in 
black solid lines, which are htted well with an exponential function (black dashed line). The 
relation of data set II is illustrated with a red dotted line. 
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165.30±15, -3.23±0.87, and 17.22±12 under 95% confidence interval, respectively. The 
trend of K of data set II is similar to the corresponding 7 . The K value decreases from 
200±31 to 69±15 km^ s“^ except the first K value is smaller than the second, and the fifth is 
also smaller than sixth. It can be seen that the K values are distinctly higher than those of 
data set I. In equation ([T]) and (E]), 7 relates to the acceleration and C relates to the initial 
velocity. The C values deduced from equation ([T]) of data set II are greater than those of 
I. It is in agreement with the previous studi es, where the horizontal velocity of GBPs is 


attenuated in AR s compared to QS regions flBerger et ah 


1998 


Mostl &: Hanslmeie 


2006 


Keys et al 


20111). This is the main reason why the K values of data set II are distinctly 


higher than those of data set I. 


4.2. DDI Approach 

Figure [H] (a) shows the relation between 7 and B in data set I using the DDI approach 
in the two-dimensional histogram. The 7 value of each GBP is measured by the slope of its 
own square displacement and lifetime on a log-log scale. The distribution of 7 in each bin is 
fitted with a Gaussian function well and the marginal distribution of 7 is projected on yoz 
plane. All distributions have similar shapes, but shift with the increasing longitudinal field 
strengths. The ranges of 7 values in different bins have no significant difference. About 
83%~90% 7 values range from 0 to 4. The mean 7 and the associated standard deviation 
are calculated by fitting the Gaussian distribution of each bin. We find that the mean 7 
values continue to mostly decrease from 1.66 to 1.33, and the associated standard deviations 
limit in a small range from 0.04 to 0.07. At the first longitudinal field strength bin, the 7 
value is 1.66±0.05, then it decreases with the increasing longitudinal field strength, and 
finally reaches a value of 1.33±0.07. On the upper xoy plane, the mean 7 values of all bins 
and a quadratic curve fit of the mean 7 values are drawn in red. 
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Fig. 6 .— The two-dimensional histograms of the diffusion indices of GBPs and the longi¬ 
tudinal magnetic held strengths by the DDI approach. The marginal distribution of the 
longitudinal magnetic held strengths is projected on the xoz plane. The marginal distribu¬ 
tion of the dihusion indices is projected on the yoz plane. On the upper xoy plane, the mean 
7 values of all bins, and a quadratic curve ht of the mean 7 values are drawn in red. Panel 
(a): the histogram of data set I. Panel (b): the histogram of data set II. 
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Figure [ 6 ] (b) shows the relation between 7 and B in data set II. The longitudinal field 
strengths range from 100 to 450 G. The distribution of 7 in each bin is also htted well with 
a Gaussian function. In different bins, about 90%~95% 7 values range from 0 to 4. The 7 
values decrease from 1.84±0.07 to 1.37±0.14, except for the value of 1.68±0.07 in the sixth 
bin. 

To explore the relation between 7 and B more clearly, we redrew the mean 7 values and 
the associated standard deviations of data set I and II using the error bars in Figure [7] (a) 
with a black solid line and red dotted line, respectively. The trend of the mean 7 values of 
data set I is htted well with an exponential function (black dashed line). By weighted curve 
htting, the empirical formula was deduced as equation ([ 2 ]), where the parameters a, h and c 
are 0.32±0.07, -2.15±1.58, and 1.41T0.06 under the 95% conhdence interval, respectively. 

Note that the marginal distributions of B of data set I and II are projected on the 
xoz plane in Figure [ 6 ] (a) and (b), respectively, which are htted with a double log-normal 
distribution with two peaks at 214±72 and 662±91 G, and a log-normal distribution with 
a peak at 277±49 G. Interestingly, there are two peaks in data set I, but only one peak in 
data set II. The hrst peak value of data set I is close to the peak value of data set II. 

The K value of each GBP is calculated by equation ([3]) using its lifetime as the 
timescale r. The mean K values are determined by the distributions of the K values of 
diherent bins, respectively. Figure [7] (b) shows both relations between K and B of data 
set I and II. The K value of data set I decreases exponentially from 89±5 to 41±4km^s“^. 
The empirical formula follows equation (|1]), where a is 66.21T11.01, b is -1.57±0.94, and c 
is 32.80±15.32 under 95% conhdence interval. The K value of data set II decreases from 
139±6 to 67±24km^s“^ except for the value of 124±9km^s“^ inside the hrst bin. We hnd 
that the K values in the hrst bins of both data sets are small, although the corresponding 
7 values are not. The main reason is that the lifetimes of GBPs with weak held strengths 
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are shorter than those with strong ones. 


5. DISCUSSION 

We nsed high spatial and temporal resolntion G-band images and simnltaneons NFI 
Stokes I and V images acqnired with Hinode /SOT. The corresponding longitndinal 
magnetic held strength of each GBP was extracted from the calibrated NFI magnetogram 
after carefnlly aligning these data with the G-band images. The point-to-point method is 
feasible to investigate the diffnsion regimes of magnetic hnx tnbes at different longitndinal 
magnetic strengths. 


5.1. Lagrangian and DDI Approach 

The Lagrangian approach and the DDI approach have been adopted to measnre the 7 
and K valnes separately. The relations between 7 and B and between K and B are both 
htted well with exponential fnnctions no matter which approach was nsed, althongh the 
valnes are somewhat different. 

The traditional Lagrangian approach is typically nsed to analyze the diffnsion of 
tracers within flnid flows. The step of mean-sqnare displacement averages the displacement 
of individnal GBPs, and then resnlts in diminished displacement. It is the main reason 
that the 7 valnes are smaller than those of the DDI approach. In addition, this approach 
inevitably depends on timescales. A short timescale generally resnlts in a small diffnsion 
index, and vice versa. Most GBPs have a short lifetime becanse the distribntion of lifetime 
follows an exponential fnnction and the mean lifetime is abont 150 s. Therefore, a few 
random long-lived GBPs will determine the diffnsion efficiency for long timescales, which 
is illnstrated in Fignre HI The timescale was determined varionsly in previons stndies. 
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Some authors took the whole mean-square displacement directly, while others cut part of 
the tail of mean-square displacement based on the percentage or the goodness-of-£t, etc. 
We calculated the mean-square displacements of all longitudinal held strength bins for the 
timescale < 300 s because this timescale is not affected by a few long-lived GBPs. 


Another approach is the DDI, where the mean diffusion index is obtained from 
the distribution of diffusion indices of individual GBPs. We analyzed all trajectories of 
individual GBPs separately. The approach avoids mixing GBPs in different diffusion 
regimes, so it sheds light on the intrinsic property of their proper motions. However, if a 
GBP moves in an erratic or circular path, the slope of linear £t on a log-log scale would be 
very small (even below zero) or very large (greater than four) with a low goodn ess-of-ht. 


Examples of such pathological cases have been given by 


Jafarzadeh et al.l (120141) . These 


cases result in meaningless diffusion indices. The percentages of meaningless diffusion 
indices of data set I and II are 13% and 8%, respectively. In detail, the very small diffusion 
indices (below zero) and very large (greater than four) account for 9.6% and 3.5% of data 
set I, 6% and 2% of data set II, respectively. Besides that, the K value of each GBP is 
calculated using its lifetime as the timescale. That means most K values are calculated 
by very short lifetimes because the lifetime distribution of GBPs follows an exponential 
function. According to equation ([3]), a shorter timescale will give rise to a smaller K value. 
This directly leads to the K values are smaller than those calculated for the timescale r < 
300 s by the Lagrangian approach. 


Above all, each of the two approaches has advantages and disadvantages respectively. 
We prefer to adopt the DDI approach to estimate the diffusion efficiency of GBPs because 
it reflects the diffusion regimes of individual GBPs. 
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5.2. Diffusion Index and Diffusion Coefficient 


The 7 and K values decrease with the increasing longitudinal magnetic held strengths, 
and they decrease exponentially of data set I. It is suggested that in the same environment, 
strong magnetic elements diffuse less than weak elements. In addition. Figure O and 
Figure [7] indicate that the 7 and K values of GBPs in strong magnetized environments are 
less than those with the same longitudinal held strengths in weak ones. The 7 values of 
data set I and II by the Lagrangian approach are 1.30±0.09 and 1.61±0.06 for the timescale 
r< 300 s, and by the DDI approach they are 1.53±0.01 and 1.79±0.01, respectively. The 
K values using the Lagrangian approach are 56±14 and 165±41km^s“^ for r < 300 s, and 
are 78±29 and 130±54km^ s“^, respectively, using the DDI approach. 


So me authors analyzed the dihusion in spe cial regions. In networ k regions of the QS or 


in ARs, 


Berger et ah 


11199811 got 7 = 1.34±0.06. 


Cadavid et ah 


(1199911 found 7 = 0.76±0.04 


for timescales s 


Lawrence et ah 


Schriiver et al 


l orter than 22 minutes and 1.10T 0.24 for time scales longer than 25 minutes. 


(1200 


indicated 7 = 1.13±0.01. 


Wang ( 1988 1 reported iF = 150km^s 


(Il996[l argued that a corrected dihusivity of OOOkm^s ^ is in a good 


agree ment with the we 


cycle (IWang fc 


Hagenaar et al. 


reported by 


Sheelev 


(19991 


per: 


brming model for the magnetic hux transport during a solar 


199411 . By tracking magnetic features in MDI magnetograms. 


Chae et al. 


bund K = 70-90km^s The smallest K value of 0.87km^s ^ was 


(1200811 from Hinode /SOT NFI magnetograms. Based on solving 


the equation of magnetic induction, they modeled the difference of the magnetic held in 
individual pixels between two frames to measure the magnetic diffusivity. The speculated 
reason for such a small K value is that their met hod aimed at individua l pixels between two 


frames with a 10 minute interval. In QS region s. 


in the rang e of 10 0-340km^s 


Yang et al. 


Chitta et al. 


Cameron et al. 


(120 llh described K lying 


( 2012 1 reported 7 = 1.59 and 77=9 0 km^ s 


( 2 OI 5 I) proposed 7 =1.50 and iF=191±20km^s 


Jafarzadeh et al. 


(l2014h 
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obtained 7 of 1.69±0.08 and i^=257±32 km^ s ^ with the Ca ll H data nsing the DDI 
approach. 


Some anthors also compared the diffnsion in different regions. 


Schriiver fc Martin 


(I 1990 I) fonnd the K valnes of 110 and 250 km^ s ^ _[n the c ore of a plage region and 


snrronnding qniet regions, respectively. 


Berger et ah 


(1l998f) determined 60±llkm^s“^ 
for network GBPs by assnming normal diffusion [K = ((A/))^/4r). However, when they 
reconstructed the velocity held using correlation tracking techniques, they presented the K 
values of 77, 176, and 28 5km^s~^ in the magnetic, network, and quiet regions, respectively. 


Abramenko et ah 


(201l|) used the TiO data obtained with the New Solar Telescope (NST) 


of the Big Bear Solar Observatory (BBSO). They indicated that 7 and K are 1.48 and 


12 k m^ s ^ for an AR, 


CH. 


Giannattasio et al 


. 53 an d 19km^s ^ for a QS region, and 1.67 and 22km^s ^ for 


(120141) dehned a network area and an internetwork area in the QS 


using Hinode /SOT NFI magnetograms. They found 7 = 1.27±0.05 and 1.08±0.11 in the 
network area (at smaller and larger scales, respectively), and 1.44±0.08 in the internetwork 
area. They also found K value ranges from 80 to 150 km^ s“^ in a network area, and from 
100 to 400 k m^ s~ ^ in an internetwork area. However, some authors had different opinions. 


Keys et al. 


(120141) estimated 7 values of ~1.2 and K values of ~120±45km^s ^ for three 


subhelds of varying magnetic hux densities. They proposed that the diffusion regimes in 
all three subhelds were the same. Their 7 and K values are consistent with those other 
authors analyzed in network regions, and more like the saturation value that we found at a 
strong longitudinal magnetic held strength. 


Previous works suggested that magnetized elements dihuse more slowly in strongly 
magnetized regions than in less magnetized ones, although the range of reported 7 and 
K values are large. Such a large range may be due to diherent instruments, data, and 
methods, especially magnetic hux densities of environments. Our results agreed with 
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their conclusions. The absence of strong magnetic helds in the medium makes it possible 
for the bright points to diffuse faster (perhaps because of less interactions), resulting in 
larger diffusion indices and diffusion coefficients. Importantly, we hnd that strong magnetic 
elements diffuse less than weak elements, both in strong magnetized environments and in 
weak ones. 


5.3. Diffusion Regime 


The magnetic field is ubiquitous in the photosphere and interacts with convective flows 
at different scales. Strongly magnetized elements are not easily perturbed by convective 
(supergranular, mesogranular, granular) motions. They withstand the perturbing action 
of convective flows much better than the wea k magnetic element s , so t heir motions are 


slower and result in smaller diffusion indices. 


Giannattasio et al.l (120141) interpreted that 


magnetic elements are transported effectively by convection, where a weak-held regime 
holds, especially internetwork; whereas, in network, where is a strong field, magnetic 


elements cannot 


De fur 


Abramenko et al. 


her transported and concentrated due to reduced convection. 


(120111 ) indicated that GBPs in strong magnetic helds are so crowded 


within narrow intergranular lanes when compared with a weak magnetic environment, 
where GBPs can move freely due to a lower population density. 

Additionally, the difference has bee n also interpreted a s the combined effect of granular. 


mesogranular, and supergranular hows. 


Spruit et al. 


decrease when approaching mesogranular boundaries. 


(Il990l) showed that me s ogran ular hows 


Orozco Suarez et al. 


(201211 found 


that magnetic elements in the internetwork start accelerating radially outward at the center 
of a supergranular cell, and decelerate as they approach the boundaries of super granules. 


Jafarzadeh et al. 


(120141) analyzed that this is caused by an increasing velocity with radial 


distance of the supergranular how prohle due to mass conservation assuming a constant 




































25 


upflow over most of supergranules. Supergranular flows systematically advect all GBPs 
toward the boundaries of supergranules, and granular and mesogranular flows impart 
the GBPs with additional velocity. The GBPs at the internal of supergranular cell tend 
to follow a super-diffusive regime because the motions of GBPs are imparted mainly by 
supergranular flow and intergranular turbulence. Once these GBPs reach the boundaries of 
supergranules, they decelerate because they are trapped in the sinks due to inflows from 
the opposite directions (i.e., from neighbouring supergranules). Granular motions are still 
active, however, causing GBPs to undergo normal or even sub-diffusion processes. 

Note that, this study only involved isolated GBPs. Non-isolated GBPs, which are most 
likely located at stagnation points, were excluded due to the difficulty of measuring of their 
displacements. Therefore, sub-diffusive GBPs might be underestimated. 


6. CONCLUSION 

We have presented a study of the dispersal of GBPs at different longitudinal magnetic 
field strengths. Two different environments are considered, namely, a strongly magnetized 
AR and a weakly magnetized quiet Sun region, characterized by different mean longitudinal 
magnetic held strengths of 132 and 64 G, respectively. The corresponding data sets were 
acquired with Hinode /SOT comprising FG G-band hltergrams (BFI) and Stokes I and V 
images (NFI). After identifying and tracking GBPs in the G-band images, we extracted 
the corresponding longitudinal magnetic held strength of each GBP from the co-aligned 
and calibrated NFI magnetograms, then categorized the GBPs into diherent groups by 
their strongest longitudinal magnetic held strengths during their lifetimes. The Lagrangian 
approach and the distribution of indices (DDI) approach were adopted separately to explore 
the dihusion efficiency of GBPs in diherent longitudinal magnetic held strength groups. 
The values of the dihusion index by the Lagrangian approach are generally smaller than 
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those by the DDI approach. The main reason is that the step of mean-square displacement 
in the Lagrangian approach averages the displacement of individual GBPs and results in 
diminished displacement. 

We hnd that the values of the diffusion index and the diffusion coefficient both decrease 
exponentially with the increasing longitudinal magnetic held strengths. The empirical 
formulas deduced from exponential htting and the parameters have been presented in 
Section 4. Stronger elements show lower diffusion indices and diffusion coefficients both in 
active regions and in QS regions. Additionally, the diffusion indices and coefficients are 
generally larger in regions of low mean magnetic hux (i.e., in the quiet Sun). 

The different diffusion regimes of GBPs are mainly set by convections. Strongly 
magnetized elements are not easily perturbed by convective motions, and vice versa. The 
reason they have slower motions and smaller diffusion indices and diffusion coefficients is 
that they withstand the perturbing action of convective hows much better than weak ones. 
Their magnetic energy is not negligible compared with the kinetic energy of the gas, and 
therefore the hows cannot perturb them so easily. 
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Fig. 7.— Panel (a): the relations between the diffnsion indices of GBPs and the longitndinal 
magnetic field strengths by the DDI approach nsing error bars. Panel (b): the relation 
between the diffnsion coefficients of GBPs and the longitndinal magnetic held strengths. 
The mean valnes and the associated standard deviations of data set I and II are drawn with 
a black solid line and red dotted line, respectively. Both relations of data set I are htted 
very well with exponential fnnctions (black dashed lines). 
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